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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-652 

THEORETICAL EVALUATION OF THE PRESSURES, FORCES, AND 
MOMENTS AT HYPERSONIC SPEEDS ACTING ON ARBITRARY 
BODIES OF REVOLUTION UNDERGOING SEPARATE AND 
COMBINED ANGLE - OF - ATTACK 
AND PITCHING MOTIONS 
By Kenneth Margolis 

SUMMARY 


Equations based on Newtonian impact theory have been derived and a 
computational procedure developed with the aid of several design- type 
charts which enable the determination of the aerodynamic forces and 
moments acting on arbitrary bodies of revolution undergoing either sepa- 
rate or combined angle- of- attack and pitching motions. Bodies with 
axially increasing and decreasing cross-sectional area distributions axe 
considered; nose shapes may be sharp, blunt, or flat faced. The analysis 
considers variations in angle of attack from -90° to 90° and allows for 
both positive and negative pitching rates of arbitrary magnitude. The 
results are also directly applicable to bodies in either separate or 
combined sideslip and yawing maneuvers. 


INTRODUCTION 


The purpose of the present paper is to develop a method based pri- 
marily on Newtonian impact theory for calculating the aerodynamic forces 
at hypersonic speeds acting on arbitrary bodies of revolution undergoing 
either separate or combined angle- of -at tack and pitching motions. Bodies 
with axially decreasing as well as axially increasing cross-sectional 
area distributions are considered in the analysis. Nose shapes may be 
sharp, blunt, or flat faced, and forebody shielding effects on flared 
aft regions are taken into account. The angle- of -at tack range considered 
is -90° 5 a ^ 90° and pitching velocities, both positive and negative, 
of arbitrary magnitude are allowed. Thus, the formulation of the problem 
is broad in scope and allows for general application of the results 
obtained to missiles and aircraft configurations of current interest for 
a wide range of flight conditions. 
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Equations are derived and a computational procedure developed with 
the use of several design-type charts which permit the calculation of 
the normal and axial forces and pitching moment for any desired body. 
Some consideration is also given the problem of modifying the results 
obtained to account for the effects of Mach number variation and differ- 
ences in the ratio of specific heats. 

Although the pitching motion referred to herein is a steady-state 
effect, it is known that the impact theory when applied to uniformly 
accelerated motion, that is, d effect, predicts a zero resulting force 
and moment. Therefore, the analysis and results of the present investi- 
gation can be considered to include the unsteady pitching condition for 
low-frequency oscillations. Because of symmetry considerations, the 
results are also directly applicable to arbitrary bodies of revolution 
undergoing either separate or combined sideslip and yawing motions, and 
low-frequency lateral oscillations. 


SYMBOIS 


x,r,o> 

R(x) 


R'(x) 


Ro 
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P 


cylindrical coordinates (see sketch (a)) 
rectangular coordinates (see sketch (d)) 
body radius in r,o> plane (see' sketch (a)) 

dR(x) 


axial rate of change of body radius, 


dx 


body radius at nose 

radius of cylindrical section (see sketch (b)) 
Mach number 

ratio of specific heats 

local pressure 

free -stream pressure 


free-stream velocity 
free- stream density 
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C N 


C L 


CD 


Cm 


G(x) 
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coefficient of pressure difference. 


P - Pn 


— p V 2 

2 P oo 00 


angle of attack 

angular pitching velocity (see sketch (a)) 
local normal component of resultant velocity 

reference area used for nondimens ionali zing purposes 
reference length used for nondimens ionali zing purposes 
body length 

Normal force 


normal-force coefficient, 


|PcoVco 2 Sr 


lift coefficient. 


Lift 


axial-force coefficient, 


ip V w 2 S r 
2 00 1 

Axial force 




drag coefficient. 


Drag 


ip V 2 S r 

2 H co 00 r 


pitching-moment coefficient, ?il°hing moment 




location of pitch axis and moment reference, measured 
positively rearward from nose (see sketch (a)) 

function defined by equation (k) 

maximum percent shielding due to forebody, defined by 
equation (25) 

maximum cross-sectional periphery at a given axial station 
shielded by forebody (see sketch (c)) 
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0 

% 


ll, 1 2 


angle subtended by periphery s (see sketch (c)) 

angle used for evaluating forebody-shielding effect, 
|o^| = i(jt - 0), (see sketch (c)) 

used as integration limits to denote portion of body 
length l 


Il[G(x)] 

l2[G(x)]l 

I 5 [G(x)] 

l[G(x)]^ 

Jl[G(x)] 

J 2 [g(x)]1 

J 5 [g(x)] 

j[G(x)]^ 


functions associated with evaluation of axial force 
(see eqs. (10) to (12) and (l6) to (l8)) 


functions associated with evaluation of normal force and 
pitching moment (see eqs. (13) to (15) and (l6) to (l8)) 


AC N " 

AC X 

AC m 

Al[G(x)] 

Aj[G(x)]^ 


increments or decrements in previously defined quantities 


parameter used in evaluating flat-faced nose contributions, 
cos a 



Superscript: 

* denotes basic contribution plus forebody-shielding effect 




k1 H OJ OJ IA 



Subscripts: 


f denotes contribution of flat-faced nose section 

t denotes total contributions from all sources 

Stability derivatives are discussed in the appendix and are defined 
therein. 


ANALYSIS 


The Newtonian impact theory is strictly applicable at Mach numbers 
approaching infinity and for values of 7, the ratio of specific heats, 
approaching unity. Studies of this approximate theory and simple modi- 
fications thereof (for example, refs. 1 and 2) have shown the Newtonian 
concept to be a valuable tool in predicting general trends and in many 
instances yielding quantitative results of acceptable accuracy when 
applied to high supersonic and hypersonic speeds. 

Basically, the theory assumes that the shock wave lies on the body 
surface and that the flow, upon striking the body, loses all of its 
momentum in the direction normal to the body surface and continues along 
the surface with its tangential component of momentum unchanged. Fur- 
thermore, only those surface areas under direct impact from the air- 
stream, that is, surface areas which "see" the flow, will experience a 
pressure force. The Newtonian approximation does not predict any effect 
on those surface areas facing away or "shielded" from the flow} in these 
regions the flow is assumed to be completely separated and the local 
pressure is equal to the free-stream pressure. 

The Newtonian hypothesis yields a local pres stir e coefficient which 
may be expressed as (see refs. 2 and 3) 



where p , p^, and are the free-stream density, pressure, and 

velocity, respectively, and where p and Vjj are the local pressure 
and local normal component of resultant velocity, respectively. It is 
to be noted that equation (l) is applicable only to those surface areas 
that face the flow; for surface areas that are shielded from the flow, 

P = P ro and C p = 0. 
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Consider the body of revolution defined with respect to the cylindri- 
cal coordinate system shown in the following sketch. 



x = Xq measured positively rearward; positive sense of rotation is 
nose-up as indicated. 

The local normal velocity Vjj resulting from a combined angle-of- 
attack and pitching motion ^ a ^ -» < q < ooj applicable to sur- 

face areas that are subject to compression flow, that is, that "see" the 
flow has been derived (refs. 3 and 4 ) and is expressible in the following 
form: 


r 


* 


V N (cn,x) 


Veo cos a 


( 


t'(x) 


1 + R ,2 (x) 


- tan a sin cd - 


q sin to 
V,,, cos a 


E 


x Q + R(x)R 



( 2 ) 


where R(x) and r'(x) define the body shape and its rate of change 
with axial distance. 


r— I CVJ CVI IT\ 
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Substitution of equation (2) into equation (l) yields the local 
pressure coefficient 

2 . I 2 

Cp(co,x) = 5 < R' (x)cos^a - R'(x)sin od sin 2a + sin 2 <n sin 2 a 

1 + R ,<£ (x) 1 

- 2 -2- sin a> cos a R' (x) Jr(x)R’ ( x) + x - x Q j 


+ 2 — sin 2 a> sin a|R(x)R'(x) + x - xjj 
Vco 

/ \2 

+ j— j sin 2 CD^R(x)R’ (x) + x - xjl \ 

\ 00 / l 


( 3 ) 


Considerable simplification is afforded by use of the substitution 


G(x) = 


-(R(x)R'(x) + x - xj) + sin a 


R'(x)cos a 

which enables equation (3) to be rewritten as follows: 




Cp(co,x) = — — - — l CQS - - a |g 2 ( x ) s in^co - 2G(x)sin o> + l] 
1 + R' 2 (x) 


or 


C p (o),x) 


gR- 2 U)cosgaj- _ G( ln n2 
1 + R' 2 (x) L J 


( 5 ) 


where G(x)sin to = 1 for R’(x) > 0 and where G(x)sin tx> ^ 1 for 
R'(x) < 0. The conditions on G(x)sin (u insure that the surface region 
under consideration is not shielded from the flow. When the surface is 
shielded, i.e., G(x)sin cs > 1 for R'(x) >0 or G(x)sin cu < 1 for 
R'(x) < 0, Cp(co,x) is, of course, zero. An examination of equation (2) 

in conjunction with the substitution of equation (4) and the requirement 
that Vu ^ 0 readily yields the indicated inequalities. It should be 
noted that, for R’(x) = 0, the parameter G(x) is infinite and a more 
convenient expression for the pressure coefficient Cp(co,x) is directly 

obtainable from equation ( 3) • 
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When the body of revolution has a flat circular face, the face 
pressure predicted by the Newtonian theory may be obtained from equa- 
tion (5) by the following substitutions therein: x = 0, R(x) = v, 

or 

R'(x) oo , and G(x) = . The formula for this pressure is given 

Voo cos a 

in coefficient form by: 




= 2 costal 1 - 


qr 


sm oo 


V, 


cos a 


( 6 ) 


and is applicable to those 

. Voo cos a 
i.e., r sm oo s . 


regions that are not shielded from the flow, 
For shielded regions, i.e., 


Voo cos a 

r sin oo > , C-p,f 


is zero. 


Another effect that requires consideration is the shielding of 
flared aft regions of the body due to the presence of the forebody. For 
example, the lateral surface area of a conical frustum would be com- 
pletely exposed to the flow at a sufficiently small a and for q = 0 
if the frustum was used as the forward section of a body; if, however, 
the same frustum was positioned as a flared aft section behind a 
cylindrical midsection, then part of the frustum would be shielded 
from the flow. A later section of the paper treats the problem of 
estimating corrections to the calculated forces and moments so as to 
account for such shielding effects. 


In the calculation of the axial force, normal force, and pitching 
moment for an arbitrary body of revolution, there are therefore three 
components to be considered: (a) the basic forces and moments associated 

with the body, exclusive of forebody shielding effects; (b) corrections 
necessitated by forebody shielding effects; and (c) the contribution of 
a flat-faced nose shape. 


The basic forces and moments (item (a) above) may be expressed in 
the following manner: 


I 


C N “ ~ 


2 

S r 


/7 

Jo J a 


0>a(x) 
o>l(x) 


R(x)Cp(to,x) sin co dm dx 


(7) 


vn ro ro m 
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pi p^(x) 

Cx = — / / R(x)R' (x)C p (ui,x)da> dx (8) 

Jq Jq>]_(x) 


pi pCDgCx) 

= — — / / R(x)Fr(x)r'(x) + x - xJlc p (o),x)sin a> du> dx (9) 

S r i r Jo Ja^Cx) u J 

where Z is the overall body length, S r and Z r are the reference 
area and reference length, respectively, used for nondimensionalizing 
purposes, and functions ojp(x) and o>p(x) are the cjo limits determined 
by the Vj^ = 0 velocity boundary and the applicable range of to, i.e., 

- ^ ts <d % Because of symmetry considerations, the multiplicative 

factor of 2 accounts for the forces and moments acting on the complete 
body. 


For body regions that are completely exposed to the flow, 
co]_(x) = and o>2(x) = for body regions where part of the lower 

surface is shielded, o>j_(x) = arc esc G(x) and u£ = and for body 

regions where pant of the upper surface is shielded, cjd^(x) = - ^ and 


^(x) = arc esc G(x). Inasmuch as the calculations of Cjj, C X> and C m 

involve (see eqs. (7) to (9)) integrations with respect to the variable 
cn of the quantities Cp(o>,x) and Cp(<n,x) sin o>, it is convenient to 
define the following functions of G(x) : 


.[GW] - / 


*/2 


arc esc G(x) 


C D (co,x)(to 


^R'(x) = 0: G(x) < -l\ ( . 
R'(x) < 0: G(x) = 1 J ^ 


sit/2 


I 2 [g(x)] = J C p ((JD,x)d£o (r'(x) >0: -1 ^ G(x) < l) (ll) 


pare esc G(x) 

(g(x)J = / ( Cp(o),x)d£D 

J -It/ ‘ 


t/2 


R'(x) ^ 0: G(x) > 1 


R'(x) < 0: G(x) = -1 


( 12 ) 


£g(x)] 


C-n(o),x)sin 0) da) 


arc esc 


R'(x) ^ 0: G(x) < -1 

R T (x) < 0: G(x) ^ 1 


J 2 [g(x 3 = f C p (co, x)sin co dco 

J -«/2 


(r’(x)> 0: -l<G(x)<l) (14) 


£(*)] = 


,arc esc G(x) 


Cp(co,x)sin co dco 


R'(x) £ 0: G(x) > 1 

R'(x) < 0: G(x) $ -1 


where the ranges of R’(x) and G(x) variations are indicated paren- 
thetically next to each integral. 

Upon integration, and designation of the various I and J functions 
as simply ljG(x7j and jJg(x)J, the following formulas and relations 

result: 

For 0 5 G(x) 5 1 and R’(x) > 0, 


- + ^ sec^a l[p(xj] = rt ( g 2 (x) + £\ 
R ,2 (x) 

1 + R ,2 (x) j |q( x j] = -2«G(x) 

R ,2 (x) 


For G(x) > 1 and R'(x) >0 or for G(x) = -1 and R'(x) < 0, 


— R,2 i x ) sec 2a I (g(x)J = [2 + G 2 (x)J 


R ,2 (x) 


arc sec 


- oW + 5 F£r aa(x) - 1 


— + ^ — — sec2a jjlj(x)] = -2G(x)arc sec-G(x) - — ^ 2 ^G 2 ^ 


x) - 1 
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For G(x) < 0 and R‘(x) >0 or for G(x) > 1 and R'(x) < 0, 

l[G(x[] = l[-G(x)) "1 
j[G(x[) = -j[-G(x]j | 

It should he noted that, for R'(x) < 0 and -1 < G(x) < 1, the 
surface region is shielded from the flow ^Cp(<n,x) = Oj and therefore 

both l[G(x)] and j[g(x)J are zero. For the special cases where 
G(x)— >±oo, a condition which arises when R'(x) = 0 and/or |a| = p 

the following formulas are applicable: 

For G(x) -> ±oo and R'(x) = 0, 


i[g(x jj = |j^-(x - x 0 ) + sin aj 2 

g[g(x)] = + ^f-p(x - Xq) + sin a 
J I *00 


For G(x)-> too, |a| = |, and R’(x) > 0 or for G(x)^+«,, 
and R' (x) < 0, 


a = 


[g(x)] = 


2|l + R ,2 (x)| l V “ 


[r(x)R*(x) 


+ X - X, 


,] + r} 


[<5(x) ^ ’ * Tj/T^wJ {t[ E(x)B ' (x) + x H + r) 


Equations (7) to (9) may be rewritten in terms of the iJ^G(x) 
and J^G(xjj functions: 


Cn = - — f R(x) j[g(x)] dx 


( 21 ) 
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C X = 


:>r o'Q 


J R(x)R'(x) I Jg(x)] dx 


( 22 ) 


J' R(x)^R(x)R'(x) + x - Xoj j[g(x)] dx (23) 


For any given body, therefore, the integrands of equations (2l) to 
(23) may readily be computed, for various x-stations and a simple graphical 
or numerical integration performed to yield the forces and moments. In 
order to facilitate such calculations, design-type charts for the I^G(x)j 

and j[G(x)] functions, computed from equations (l6) to (l8), are pre- 
sented in figures 1 and 2, respectively, which are contained in a pocket 
attached to the hack cover of this report. 

The second component to he considered in the force and moment cal- 
culations is (item (h)) the correction necessitated hy forehody shielding 
effects if applicable to the given body for the specified angle of attack 
and pitching velocity. In this connection, it should he noted that an 
exact correction would he a function of many variables and thus would 
not he useful as a simple design- type procedure. It is also felt that 
for most bodies this particular effect is not too significant overall and 
thus a calculation to define the maximum correction possible for a simple 
shape would he more appropriate. 


Consider, then, a flared aft section situated behind an infinite 
circular cylinder of radius R^ shown in the following sketch: 


I 




Sketch (c) 


Sketch (h) 
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At an axial location x behind the cylinder, the maximum cross-sectional 
periphery shielded by the cylinder is s (see sketch (c)) where 


s 


0R(x) 


R(x)arc 


cos 



or by trigonometric substitution. 


s = 2R(x)arc cos 




(24) 


The quantity 


expresses the maximum percent shielding that may 


2itR(x) 

be present due to the forebody. The following relation exists: 


s , 100 1 

— - = k = arc cos | 

2rtR(x) « 


Ry 


V R ( X )/ 


(25) 


and it is readily apparent that the maximum percent shielding ranges 


/ K b \ / R b 

from 0 = 0 to 50 percent 

\R(x) J \R(x) 


= 1 ) . The variation of the param- 


eter k with the ratio 


R> 


R(x) 


is shown in figure 5* 


The parameter k may in turn be related to the integration limit 


k = 


2R(x)(| - a*) 


2itR(x) 


X 100 = 50 - ^ cot, 


(26) 


where 


cd^ = arc sin 111 - 


' R b ; 

V R ( X )> 


(27) 



Omitting calculation details, the procedure is to derive expressions 

for Al[G(x)j and Aj[G(x)j through use of equations (10), (12), (13 ) , 

and (15) with appropriate changes in the 00 limits of integration so as 
to account for the shielding boundary Corresponding corrections to 


the force and moment are then introduced: 


C N* = C N + ac N 


C X “ °X + 


Cm = C m + AC, 


A CX = ‘ ^ / Z 2 R(x)R ' (x) M [ G H ** 


AGnt = - 


r ^ Aj[g(x)] dx 


ACm = — y" J 2 R(x)[r(x)R'(x) + x - x Q Aj[g(x)] 


In equations ( 29 ) to (31 ) > the integration is performed over that sec- 
tion of body length affected by the forebody, herein denoted from Zq 


to Z 2 . The functions AI 
follows : 


[g(x>] 


and AJ 


[GW 


are tabulated as 



vji ro ro h tr 1 
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For 0 = G(x) ^ 1 and 0 ^ cub = p 


— + — — — — sec^a Ai[g(x) 1 = G^(x)^ - aj^ + cos sin 

R , 2 (x) V ' 


~\ 


1 + R * l? ! sec2a Aj[g(x)] = -2 
R' 2 (x) 


- 4G(x)cOS 0)Jj + Jt - 203b 
G 2 (x)cos 03^(2 + sin^03b) 


G(x)^03b - cos 03|j sin % ” + cos % 




(52) 


For 1 ^ G(x) < 00 and 0 ^ arc esc G(x) 

^ - - ^ sec2a AI Jg(x)J = |g 2 (x) + 2jaxc sec-G(x) 


"\ 


1 + R ,2 (x) 


+ 3\/g 2 (x) - 1 - 4G(x)cOS 03b 
+ G 2 (x)^COS 03b sin 03b - “ - co^ - 203b - Jt 

— sec2a AjJg(x)J = -2^cos a3b + G(x)^ + 03b - cos 03b sin co^) 

G 2 (x)cos 03b (2 + sin 2 ab) 

3 

- G(x)arc sec-G(x)j 




( 33 ) 
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For -1 % G(x) < 0 and ^ 0, or -oo < G(x) < -1 and 

arc esc G(x) < S 0, 


Al[-G(x)] =Al[G(x)] 

Aj[-G(x)] = -Aj[g(x)] 


For any other combinations of the parameters and G(x), the fGre- 

body does not give rise to any additional shielding and the correction 
factors are zero, that is, AI^G(x)j = 0 and AJ^G(x)j = 0. The varia- 
tions of AI [g(x)] and AJ^G(x)j with the parameter G(x) are shown 
in figures 4 and 5 for several values of maximum percent shielding. 
Included in these figures are the basic lj^G(x)j and J^G(x)j curves 

so that the relative importance of the corrections can be assessed. The 
corrections AGr, ACx, and ACm may now be readily obtained for any 

given body in accordance with equations (29) to (31) by using the identi- 
cal procedure discussed previously for calculating Cjj, C^, and C m . 

Although the corrections AI [G(x)j and AJ {g(x)J just discussed 
are directly applicable to flared aft sections situated behind forebodies 
for which R'(x) ^ 0, a similar approach utilizing an average or modified 
value of R^ could readily be applied (if required) to body shapes with 

regions characterized by R'(x) < 0. 

In the event the body of revolution considered has a flat face, a 
third contribution (item (e)) to the force and moments is involved. The 
incremental normal force will, of course, be zero but the axial-force 
and pitching -moment increments will be finite. 

The pressure acting on the flat face is given by equation (6) as 


sin a) 


cos a 



Consider the following sketch: 



If, for convenience, the parameter a is introduced to denote the 

COS (X 

quantity — - — , it is readily seen that no shielding effects will be 


present for values of | a| ^ 1 and the following integral expressions 
may be written: 



(35) 


Cm,f 


4 cos^a 
S r Z r 





2 

) z dy dz 


(36) 


For values of |a| < 1, the equations above apply provided the following 

changes in limits are incorporated to account for the shielding boundary: 
for q > 0, replace upper limit on z by aR 0 , and for q < 0, replace 
lower limit on z by aR Q . 


l8 


The indicated operations have been carried out and the following 
closed-form formulas result: 

For | a | £ l. 


c N,f = 

0 


ft — 

2 p 

itR 0 cos c a ]_ 

+ k, 

^X,f “ 

S r 

2a 2 

r ^ = 

X p 

aRq^cos a 

1 

^m,f 

s r z r 

a 


(37) 


For Ja| < 1, 

"N 

c N,f = 0 


n R 0 ^cos 2 a 3(^&2 + l)arc cos- 

|a| + | a J ( 2a 2 + 13 )|fT - a 2 

^’ f s r 

6a^ / 

„ _ E 0 Jcos2a 15aarC COS -| a ! 

" -ffr( 2a ■ 9a2 ■ ^i 1 ■ a2 

m,f s r z r 

15 a 2 J 


The variations of C x ^-f and f with the parameter a are 

shown in figure 6. It should be noted that only positive values of a 
need be considered since C x ^ f (a) = C x ^f(-a), and C^fCa) = -C m ,f(-a). 

Thus the three contributions to the forces and moments have been 
considered: (a) the basic components Cy, (b) the forebody 

shielding components ACjj, AC X , AC m , and (c) the flat-face contri- 
butions C m ^. These components are additive, 

°N,t = C N + ACjj + Cj^f 
C X,t = C X + ^x + C X,f\ 

^m,t = C m + ACm + C m ^f 

J 


s 


(39) 


hi rH OJ CVJ LT\ 



t-1 f-H OJ CVI tf\ 
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In the event the lift coefficient C L and drag coefficient Cp 
are desired, they may readily be obtained by the formulas: 


C L = C N cos a " C X sin a ' 
Cd = Cjj sin a + Cx cos a 


(^ 0 ) 


For values of a and q in the immediate vicinity of zero, a con- 
venient form for evaluating stability derivatives may be obtained from 
the equations derived herein; they are presented for general interest 
in the appendix. 


APPLICATION OF METHOD 


In order to illustrate the procedures used in evaluating the aero- 
dynamic forces and moments, consider the body of revolution shown in 
figure 7 consisting of a hemispherical nose faired into a cylindrical 
midsection followed by a flared aft region composed of two conical 
frustums. The equations defining body shape and slope are given in 
the figure along with values chosen for aerodynamic parameters ^ angle 

of attack a = 20°, positive pitching — = 0.015 about center of 


x 4\ 

gravity -2. = —J, For purposes of convenience, the body length is 

chosen as 9 units and the corresponding lengths of other quantities are 
indicated. 


Sufficient axial stations are chosen (21, indicated by the dots 
along the body axis of symmetry) and a value of G(x) computed from 
equation (4) for each value of x. For each G(x) value so obtained, 
there are read off from figures 1 and 2, respectively, values of the 
IjG(x)j and j|g(x)J parameters. (For values of G(x) -> co, simple 

computations are made directly from the formulas included in the figures.) 
Then in accordance with equations (21), (22), and (23) the functions 


R(x) •<j[g(x)] , R(x)-R'(x)-i[g(x)] , and R(x) • [r(x)R* (x) + x - x 0 ]-j[g(x) 


are tabulated and plotted against the axial coordinate x as shown in 
figure 8. Integration by any suitable means (Simpson's Rule, plani- 
meter, etc.) to obtain the areas under the solid-line curves of 


v 
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figure 8 and multiplication by the appropriate factors of equations (21), 
(22), and (23) yield the following results for the basic forces and 
moments : 


C X = 0.7^ 
C N = 0.5^ \ 

C m = -0.33 


(41) 


To calculate the correction due to forebody shielding, that is, 
AC X , AC^, and AC^, values of the parameter k defined by equa- 
tion ( 25 ) are obtained from figure 3 for those axial stations behind 

the cylindrical section. The parameters AI^G(x)j and AJ^G(x)j are 

then estimated from figures (4) and ( 5 ), respectively, and the calcula- 
tion of the incremental forces and moment defined by equations ( 29 ), 
(30), and ( 31 ) proceeds in a manner directly analogous to the procedure 
described above for the basic contributions. The incremental effects 
are shown as dashed lines in figure 8 and result in the following 
corrections: 


AC X = -O.O3"' 

ac n = 0.03 > 
AC m = -0.02 


(42) 


Inasmuch as the body chosen for illustrative calculations does not 
have a flat surface at the nose, the flat-face contributions (eqs. (37) 
and (58)) are, of course, zero. Thus in accordance with equations (28) 
and (39), the total coefficients are: 


c X,t " °*71 "] 


C N ,t = 0.57 } 


Cm,t = -0.35J 


( 43 ) 
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It should be noted that the Newtonian theory does not account for 
either variations in Mach number or for differences in the value of 7 , 
the ratio of specific heats. The theory is strictly valid at Mach 
numbers approaching infinity and for 7=1. Simple modifications to 
the Newtonian theory have been suggested to account for these variations. 
As an example for blunt bodies, the Newtonian coefficient might be multi- 


plied by the ratio (see ref. 5 ) j 


7 + 3 

2(y + 1) 


1 


2 

(7 + 3)M 2 ' 


For pointed 


bodies with attached shocks, the Newtonian values appear to be more 
satisfactory. Some discussion pertinent to Newtonian theory modifica- 
tions may be found in reference 6 . 


CONCLUDING REMARKS 


Equations based on Newtonian impact theory have been derived and 
a computational procedure developed which enables the determination of 
the aerodynamic forces and moments acting on arbitrary bodies of revolu- 
tion undergoing either separate or combined angle- of- attack and pitching 
motions. The analysis considers variations in angle of attack from -90° 
to 9 O 0 and allows for both positive and negative pitching velocities of 
arbitrary magnitude. The results are also directly applicable to bodies 
in either separate or combined sideslip and yawing maneuvers. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., October 28, i960. 
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APPENDIX 

CALCULATION OF STABILITY DERIVATIVES 


Although the Newtonian theory predicts forces and moments that 
are nonlinear with respect to variations in angle of attack and/or 
pitching velocity, considerable attention has been given the problem 
of calculating stability derivatives, that is, rates of change of forces 
and moments expressed in nondimens ional form (for example, refs. 3 
and 7). Inasmuch as the derivatives are conventionally evaluated as 
a -» 0 and q -» 0, it is apparent that considerable caution must be 
exercised in the application of results so obtained without regard to 
nonlinear effects. 

The equations for the basic forces and moments derived in the pres- 
ent paper may be readily expressed in the form of stability derivatives 
by appropriate differentiation of equations (7), (8), and (9) with the 
limits o>l(x) and ^(x) therein replaced by -n/2 and rt/2, respec- 
tively, and utilization of the evaluations for I2 |G(x)j and J2 J^G(x)j 
as given by equation (l6). The following formulas result: 


-N, 


a 


/ dcA _ ij*t C l R(x)R' 

/ a -*0 S r ^0 1 + R* 


(x) 

2 (x) 


dx 


(Al) 


q— »0 


dC- 


c N r 


N 


2V 

^ Voo / q -*0 
a ->0 


_ 8 n r l R(x)R , (x)[r(x)R'(x) + x - Xq] 
S r z r Jo 3 . + R’ 2 (x) 


dx (A2) 


■'“a 


/^Cjj^ 

| _ ^ 

pi R(x)R'(x) 

[r(x)R'(x) + x - xj 

\8a ) 

a-»0 S r Z r ^ 

q ~>0 

0 1 + R ,2 (x) 


dx (A3) 
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-mq 


„ j 

07r 

2V„ 


8rt r l R ( X ) R '( X ) [r(x)R' (x) + X - Xoj 

J n 


S r z r ^ 'JQ 


q-> 0 
CL — > 0 


1 + R ,2 (x) 


dx (a4) 




^CjA 

(da/a-*0 

q-»0 


= 0 


(A5) 



a-»0 


(A 6) 


At the limit a = 0 and/or q = 0, the normal -force coefficient Cjj 
and the pitching-moment coefficient are both zero, but the axial- 

force coefficient is finite. A formula for Cx may be derived 

through use of equations (8) and (l6) and is given by 



it* f l R(x)R ,? (x) ^ 
Sr 'JO i + R’ 2 (x) 


(AT) 


Equations (Al) to (A7) are valid for arbitrary bodies of revolu- 
tion and may be readily evaluated by analytical, numerical, or graph- 
ical procedures. 

In the event the body under consideration has a flat face at the 
nose, there are some additional contributions to be considered; these 
may be obtained by straightforward differentiations of equations (37) 
and are as follows: 


h 
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GO — > 0 


where R 0 is the radius at the nose. 

At the l imi t a = 0 and/or q = 0, the normal-force coefficient 
C N ^ and the pitching -moment coefficient C m ^f are both zero, but the 

axial-force coefficient is finite and is obtained from equation (37) as 


\ji ro ro f— 1 t -1 
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(Mo*> 

q=0 


s r 


(Al4) 


Additional increments due to forebody shielding have not been con- 
sidered) these effects are felt to be relatively insignificant for 
values of a and q at or in the immediate vicinity of zero. The 
effects of centrifugal forces in the flow have also been excluded in 
the analysis. 

In the event stability derivatives are desired for values of a 
and/or q significantly different from zero, it is suggested that force 
and moment variations with a and q be obtained by the method of the 
present paper so as to emphasize any nonlinear characteristics. The 
required slopes may then be estimated by the usual procedures. 
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Figures 1 and 2 are contained in a pocket attached to the hack 
cover of this report. 



R(x) 

Figure 3.- Variation of the parameter k with the ratio R^/p^x) 
where k, the maximum percent shielding due to forebody, is 
defined by equation ( 25 ) . 
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Figure 6.- Chart for evaluation of the force and moment contributions and 








